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APPLICATION OF THEODORSEN’S THEORY TO PROPELLER DESIGN 


By John L. Ceiglee 


SUMMARY 

A theoretical analysis is presented for obtaining by use of 
Theodorsen’s propeller theory the load distribution along a 
propeller radius to give the optimum, propeller efficiency for 
any design condition. The efficiencies realized by designing 
■for the optimum load distribution are given in graphs, and 
the optimum efficiency for any design condition may be read 
directly from the graph without any laborious calculations . 
Examples are included to illustrate the method of obtaining the 
optimum load distributions for both single-rotating and dual- 
rotating propellers. 

INTRODUCTION 

Recent contributions to the theory of propellers have been 
made by Theodorsen in a series of reports (references 1 to 4). 
In the first report of the series (reference 1) a method based 
on electrical analogy was devised for obtaining the ideal 
circulation functions for single-rotating propellers. These 
circulation functions were shown to be in good agreement 
with the theoretical calculations made by Goldstein in 
reference 5 for two- and four-blade single-rotating propellers 
and with the extrapolations to other numbers of blades 
made by Lock and Yeatman in reference 6. The electrical- 
analogy method of measuring these functions was also 
applied to more difficult cases for which no theoretical 
calculations had previously been made; in particular, to the 
case of dual-rotating propellers. 

Theodorsen in reference 1 introduced the concept of the 
mass coefficient k, which is an integrated value of the circula- 
tion functions. The mass coefficient represents the effec- 
tive cross section of the column of the medium pushed by 
the propeller divided by the projected-propeller-wake area. 

This mass coefficient is made use of in. the development of 
Theodorsen’s theory. In reference 4, expressions are given 
for computing the thrust, the energy loss, and the efficiency 
of any propeller with ideal circulation distribution based on 
the conditions in the final wake in terms of the mass coeffi- 
cient. It is of interest to mention that the mass coefficient 
or mass of air operated on by the dual-rotating propeller is 
much greater than that affected by the single-rotating 
propeller for the same set of operating conditions. This 
large difference in the mass coefficients for the two cases 
indicates that calculations for dual-rotating propellers 
based on the ideal circulation functions for single-rotating 
propellers are inadequate. 

Theodorsen’s theory, as previously mentioned, is based on 
the conditions in the final wake. The present analysis 


attempts to interrelate the conditions in the final wake to the 
propeller and to give the information necessary to design a 
propeller for any desired operating condition. For single- 
rotating propellers, the method yields the same results as the 
conventional vortex theory with the Goldstein tip corrections 
applied. By the conventional vortex theory, however, it is 
necessary to determine the optimum blade-load distribution 
and then to make element strip-theory calculations in order 
to obtain the optimum efficiency' for a given design condition. 
This procedure has been followed in reference 7 for a wide 
range of operating conditions. By Theodorsen’s theory the 
optimum efficiency ij can be obtained directly for any design 
condition from its relationship to the mass coefficient without 
laborious calculations. Thus, in the selection of a propeller 
for any design condition, a close estimate of the efficiency 
can be obtained before the design is made. 

The circulation functions and mass coefficients for the dual- 
rotating propeller were obtained in reference 1 for the ideal 
case and refer to conditions in the ultimate wake. Both 
propellers were assumed to operate in the same plane. Ob- 
viously, this condition is unattainable in the design of an 
actual propeller. The degree to which the ideal case can be 
realized in practice, or the applicability of the ideal functions 
to a given case, require further consideration and confirmation. 

SYMBOLS 

B number of propeller blades 

b chord of propeller-blade element 

c t section drag coefficient 

Ci section lift coefficient 

P c Ideal power coefficient (e,+e) 

Pc T total-power coefficient (P c -\-t r ) 

c, thrust coefficient 

c, T net thrust coefficient (c,— f tt ) 

D diameter of propeller 

d drag of propeller section 

D a diameter of wake helix surface 

E ideal energy loss in wake ^ pFiw* w-pr; V 

E d energy loss due to blade drag 

e induced energy loss coefficient 




S3 
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projected area of helix (at infinity) 
circulation function 
lift of propeller section 

propeller rotational speed, revolutions per second 

input power to propeller 

tip radius . 

radius to any blade element 

thrust of propeller 

power loss due to drag (nondimensional) 
axial power loss due to drag 
rotational power loss due to drag 
forward axial velocity of propeller 
axial interference velocity (at propeller) 
average axial interference velocity (at propeller) 
resultant interference velocity (at propeller) 
rotational interference velocity (at propeller) 
average rotational interference velocity (behind each 
propeller) 

resultant velocity on the propeller at radius r 
local self-interference velocity 

rearward displacement velocity of helical vortex 
surface 

ratio of displacement velocity to forward velocity 
(w/V) 

radial location of blade element ( r!R ) 
angle of attack, degrees 
induced angle of attack, degrees 
blade angle, degrees 

advance ratio 

geometric advance ratio (F/xrcD) 
mass coefficient ^2 j K(x)x dx'j 
axial energy loss factor 


propeller efficiency ^ 


jj( ideal propeller efficiency ( c,jP 0 ) 

p mass density of air 

a propeller element solidity (J?&/2nr) 

ae l propeller element load coefficient 

r circulation at radius a: ^ r (x) = 2 T W) VJ K{x)^j 

<j> angle of resultant velocity W at plane of rotation 

, * -i VlnD 

0o= tan 1 

oj angular velocity ._ ... 

Subscripts: 

F front 

R rear 

0.7 R at 0.7 radius 


OPTIMUM PROPELLER DESIGN 
SINGLE-ROTATING PROPELLERS 

Velocity diagram. — The velocity diagram for the single- 
rotating propeller is shown in figure 1. This figure is a 
reproduction of figure 13 (reference 2) .with some additional 
designations. The relationship between the axial inter- 
ference velocity at the radius_r, as given by the vortex theory, 



Figure l— V elocity diagram for single-rotating propeller. 

to the displacement velocity w of the vortex sheet is calcu- 
lated in reference 2 and is shown in figure 1. The forward 
axial velocity of the propeller is F and the tangential velocity 
with respect to the air at rest is or. The vector BD is the 
resultant interference velocity F< of the air with respect to 
the air at rest. Thus, the resultant velocity TF of a point on 
the propeller at the radius r is given by the vector CD. The 
lift force l is perpendicular to this vector and the drag force 
d is exactly opposite in direction to IT as indicated. From 
this figure a comparison of the method of analysis presented 
herein may be made with tho conventional vortex-theory 
method. It is required to find the point D in order to locate 
the end of the velocity vector IF and the angle 0 that the 
vector JPF makes with the direction of rotation. By the 
conventional vortex theory, the point D is located by starting 
with point B obtained from tho V/nD of the undisturbed 
flow, proceeding in the F-direction the distance F«, and 
then taking the perpendicular to this direction a distance 
V r . (See reference 8.) The angle <j> is given by 
Y+V a V+ V 

tan0 =— — XT and TF=- . - , • In the calculation of interfer- 
ur— V T sm 0 

ence velocities F tt and V r the local tip correction or Goldstein 
factor must be used to obtain the correct location of tho 
point D» 

With the method developed in references 1 to 4, only tho 

value of which remains constant with radius, need be 

used. With this concept it is possible to use the integrated 
values of the mass coefficient as determined by the electrical 
analogy of reference 1 to obtain the detailed information 
needed at any radius. By this method the point D can be 

located by proceeding from point B a distance in the 

F-direction to the point E and then down the direction of 
the velocity vector IT 7 a distance DE, where DE is obtained 

from the geometry of the figure as g w s * n 4> and 

V+±w 

tan 0= - (1) 



The resultant velocity is 
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I Using equations (3) and (5) for r gives at once the identity 


w =-^r~2 wshi<l> 


, (V J r w)w Tr , . 2 sin <f> 

be ‘ = Bn R{x) i “ 


F+jW cosV 


= ^h( v+ ^ wcos^,^, ) (2) 

The interference velocities may be obtained from the geom- 
etry of the figure by 

U<=^ w cos <t> 


and 


V a =Vi cos w cos 3 <f> 


T T r =T\ sin 4 >~ 7 j w sin <j> cos <£ 


Optimum blade-load distribution. — The design problem of 
an optimum propeller consists essentially in obta ining the 
value of the element load coefficient bc t at each radius of the 
propeller blade. With the direction and magnitude of the 
relative velocity given at each station there remains only the 
choice of a section to give efficiently such a lift at the appro- 
priate angle of attack. The value of e t should be at or near 
the ideal lift coefficient for the section in order to give 
mini mum drag coefficient. 

The method developed in references 1 to 4 treats the 
velocity w as an independent parameter upon which all the 
other quantities depend. This reversal of procedure is con- 
venient since all quantities are actually functions of w. The 
velocity w is related to the power coefficient P c of the pro- 
peller and also to the element load coefficient <rc t . The rela- 
tion of w to crc t is developed herein and the relation of w to 
P c , which must be obtained in order to use it for design, is 
given subsequently in the section “Procedure for Propeller 
Design.” 

The required ideal circulation T(x) is given in reference 1 by 
V(x)= MV ^ w)w K{x) 


(U+w)w 

Bn 


K{x) 


(3) 


In order to determine the element load coefficient bc t the rela- 
tion for the equality of the force on a vortex element and on 
an element of a lifting surface is given as 

P nr=4 pW'-cj, 

where b is the chord of the element. Hence, 

r=4 Wcfi (4) 

where IF is given in equation (2), and thus 

r= 4 ibh ( r+ ^ w cosV ) Cib (5) 


Introducing the nondimensional velocity v}=^< the 


Bb 


V +2 W 


solidity <j=— j and tan 4>= 2 (equation (1)) gives the 


nondimensional relation 

l+w 


crC|= 


(l +| w) -f-| V) cos’*) 


2 wK (r) 


sin 
cos <j> 


( 6 ) 


The selection of a propeller for a given airplane installation 
may be based on a method of evaluating a series of propellers 
for various operating conditions in order to determine the 
most suitable propeller. It is probable that several propellers, 
varying in diameter, blade number, propeller operational 
speed, and direction of rotation are equally as efficient for 
the design condition so that other considerations may enter 
into the propeller selection. However, the optimum effi- 
ciency for the propeller selected may be obtained from the 
charts, and therefore the load distribution along the radius 
that will give this optimum efficiency remains to be 
determined. 

The value of <rc t may be calculated for any radius from the 
relation 


<rCi= 


1+W ' 


where 




, . i v 1+ !” , , , I +!"’ 

0=tan 1 =tan 1 X, 

v nD x 1 x 


(7) 


DUAL-ROTATING PROPELLERS 

In the design of dual-rotating propellers, it has been 
customary to select two propellers designed for single 
rotation and to use them as a dual-rotating propeller. The 
fact that the circulation functions and the mass coefficients 
obtained by the electrical-analogy method (reference 1) 
are very much larger for the dual-rotating propeller than the 
sum of the values for the two single-rotating propellers 
indicates that the functions as used heretofore are not proper. 
The electrical-analogy method represents the case of an 
idealized dual-rotating propeller in which the two components 
are in the same plane with the same load distribution on each 
component and with equal power absorption. Since actual 
propellers cannot conform to this ideal case, the applicability 
of the ideal functions requires further confirmation. Never- 
theless, the optimum distribution for the dual-rotating 
propeller is essentially different from the single-rotating 
propeller, and in this analysis the loading functions and the 
mass coefficients as determined by the electrical-analogy 
method are assumed to apply to the optimum dual-rotating 
propeller. 
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Interference velocities for dual-rotating propellers. — The 
average axial interference velocity far behind the propeller 
obtained from the momentum considerations is 

2 V a =KW 

where k is the mass coefficient and w is the axial displacement 
velocity. This mean value is equally due to each of the two 
oppositely rotating propellers. The average axial interfer- 
ence velocity due to each is therefore exactly 

Tt 1 

Va — 2 KW 

The average interference velocity at the propeller plane is 
one-half the value in the final wake and, therefore, 



where ^ V a represents the average axial interference velocity 

at the propeller plane due to each component of the dual- 
rotating propeller. With the two propellers separated by 
a small axial distance, this velocity refers to a plane between 
the two propellers. The interference velocity at the front 
propeller is smaller and at the rear propeller is larger than at 
the plane between the propellers. In the following treat- 
ment, the propellers are considered to be very close together 
so that the axial interference velocity is the same on both 
propellers. 

In the final wake, the mean value of the rotational inter- 
ference velocity for the ideal case is given by 

2V r =0 

For an infinite number of right and left blades equally loaded, 
rotational components would cancel exactly. However, the 
average rotational interference velocity immediately behind 
each propeller may be considered as 

F r =^ kw tan 4> 


In summary, the mean interference velocities acting on the 
front propeller from the rear propeller are 
Axial: 



Rotational: 


y r =o 


The mean interference velocities acting on the rear propeller 

from the front propeller are 

Axial: 



Rotational: 


V r =^ kw tan <j> 


It is useful to recognize that the mean self-interference of 
each propeller in its own piano is 

Axial: 

1 

T KW 

4 

Rotational: 

kw tan <f> 

Velocity diagram for the dual-rotating propellers. The 
velocity diagram for the dual-rotating propellers is shown 
in figure 2. As in the case for the single-rotating propeller, 
the axial displacement velocity at the propeller is equal to 

^ w. In figure 2 the vector AB gives the mean axial inter- 
ference velocity ^ kw of each propeller acting on the other 
propeller. The vector BC gives the mean rotational inter- 
ference velocity kw tan <j> of the front propeller acting on 

the rear propeller. The total interference velocity acting 
on the front propeller from the rear propeller is therefore 
given by AB, and the total interference velocity acting on tile 
rear propeller from the front propeller is equal to the vector 
AC. The local self-interference velocity of the front pro- 
peller is given by RV, and the corresponding helix angle is 
given by fo-. The local self-interference velocity of the rear 
propeller is given by and the corresponding helix 

angle is given by 4> R . The angle <f> F is slightly larger than 
the ideal helix angle <f> given by the displacement velocity 

w and 4>r is slightly smaller than <j>. The design condition 

of most interest is the one for which I> for each blade of the 
front propeller is equal to r K for each blade of the rear 
propeller. The number of blades on the front and rear 
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propeller are considered equal and the rotational speeds, the 
same. This condition gives the self-interference velocity 
on the front propeller equal to the self-interference velocity 
on the rear propeller and means that D and E must be at 
the same horizontal level. 

As <f>r and <f>& are needed in the design of the propeller, it 
is seen from figure 2 that the associated displacement velocity 
on the front and rear propellers has been increased and 
decreased, respectively, by the amount 

Aic=i kw tan 2 * 

The displacement velocity is therefore 
Front: 


Rear: 


tan 1 *) 
\ tan 1 *) 


From figure 2, the velocity TTj- is shown to be given by the 
relationship 

TTV= ^ + I KWsin ^ 


^1 + j kw sin ! * 0 ) 


sin * 0 ' 

and the angle 4> r is given by 

w (l K tan 2 *) 


( 8 ) 


tan 4 > f ~- 


c or 




(9) 


where * is given by the relationship 
tan *= 


T— ■ 1 

F-hjw 


«r 


= WVx( l +\ w ) 


Similarly, 


V 1 1 

^ aw sin * 0+2 ^ an 0o cos 0o 

V ,3 

= -j - , +-rKW sm *0 
sm * 0 4 


and 


= sin*i( 1 "*'I m sinZ 0o) 

T 7 +-| w(l~^K tan 2 *') 
tan **= — ^ 

w r 


( 10 ) 


( 11 ) 


Optimum blade-load distribution, — The optimum distri- 
bution of blade loading is obtained from the determination of 
the element load coefficient be { at each radius from the funda- 
mental relation 

^ p&cjir 2 =prTr 

where r has been given in equation (3) by 

t "» (V+w)w T „_, 

T ~ ~Bn K{x) 

Eliminating r gives 

\ Bbc,W= (V ^ - w)w K(x) 

A 71 

but ^ ^~ r ~ (T 1® H 16 solidity of each component of the dual- 

rotating propeller, if the number of blades in each compo- 
nent are assumed to be equal. Therefore, 

<rc t W=- ( 1 +w)wVK(x) 

For the front propeller, this equation may be solved by use 
of equation (8) 


r r \ _ Z_ J_ (l+w)^sin* 0 

,1 _ . . K{x) 

1 kw sm 2 * 0 

and for the rear propeller by use of equation (10) 

r V 1 (l+w)wsin* 0 rrro 
{aC ^~W - - 3 _ . K(x) 

1-t^KW sm z * 0 


USE OF DESIGN FORMULAS 


( 12 ) 


(13) 


In order to use the relation for c re„ note that it contains not 
only the independent variable w but also the function K(x) 
and the angle *. The parameter K(x) should be expressed 
V~\~w 

as a function of — r— > which is based on the wake helix 

71L/q 

diameter. As was shown in reference 3, however, D 0 differs 
only slightly from the propeller diameter D and in the present 
design procedure D is used instead of D a . The function 

V ' w 

K{x) for single-rotating propellers is plotted against in 

figure 3. Similar plots for dual-rotating propellers were 
taken from reference 1 and are presented in figure 4. 

EQUATIONS FOR PERFORMANCE CALCULATIONS 
SINGLE-ROTATING PROPELLERS 

In reference 4 the thrust has been given by 

T= P F KW [F+io(J-j-f)] 

and the ideal energy loss in the wake has been given by 
E=pFnw t (jj w+l y) 
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(a) Two-blade single-rotating propeller (reference 6). (b) Three-blade single-rotating protieller (reforeneo 0). 


Figure 3.— Circulation function K(x) for single-rotating propellers. 
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-5 IX) /£ 2X) 3X) 


V+w 
nD 

tc) Four-blade single-rotating propeller. (d) Six-blade single-rotating propeller. 

Figure 3.— Continued. 
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0 .5 1.0 1.5 2.0 2.5 3.0 


V+w 

nD 

(6) Eight-blade aingle-rotating propeller. 

Fkhjre 3.— Concluded. 

. . . . in 

With the introduction of the nondimensional quantity w=y> 

the thrust coefficient in nondimensional form is 

T 

1 

\t>\ n F 

—2kW £l +w (l+f)] (14) 

and the induced loss coefficient is 

E 

«-T - 

| pVF 

=2«ii?(4+4w) (15) 



Figure 4.— Circulation function £T(r) for dual-rotating propeller* (reference 1). 


The power coefficient P c =c,+e is given by 

P e =2Kw(l+w)(^l-j^w'^ (10) 

The efficiency is given by 



These formulas axe all that are necessary for single- 
rotating propellers. The performance of the dual-rotating 
propeller is computed by the same formulas. 

DUAL-ROTATING PROPELLERS 

The thrust of the front propeller is given by 

dTf=4 p(2 7rr)HV (ac t ) F cos 4>f dr 

and with (<rc t ), from equation (12) and 11? from equation (8) 
T,=„ § V»(I+ffi) l £'(l +1 m sin*0o) ^ md* 

(18) 
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0 12 3 4 5 6 7 


P+io 

71 J) 

(b) Eight-blade dual-rotating propriler. 

Figure 4.— Continued. 

Similarly, for the rear propeUer 

r*=e f *"55(1+55) i [‘O+I * “■*•) ffiTW K 

(19) 

The coefficients c„ e, and P e for the dual-rotating propellers 
are given in the same form as in equations (14), (15), and 
(16) for single-rotating propellers. The only difference in 
the coefficients results from differences in the values of 
k, w, and e/ie which are substituted in the equations. 

BLADE-DBAG LOSSES 

The frietional loss or loss in efficiency due to the profile 
drag of the blade is 

E d =B | J* bCiW* dr 

The drag force per unit length is ^ pW 2 bCi where IF, the 
resultant velocity of the blade element, has been given in 



V+1D 

nfi 


(c) Twelve-blade dual-rotating propeller. 
Figure 4.— Concluded. 


equation (2) for single-rotating propellers by 


W= tih>( v+ ^ wc0 ^) 


For the design condition, w is small, and because of the 
obvious uncertainties in the determination of the value of c it 

it is not necessary to retain the second term ^ w cos 1 <j>. In- 
troducing the solidity factor tr=^^, permits the drag loss to 
be given by 


E d =tR* p V* P -^hxdx 

Jo sm a 0 


or in nondimensional form 
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The component power losses are then, to the same degree 
of approximation in nondimensional form, 

Rotational: 


f = 2 f 1 a* d x 

T X/Jo sin 4> 

(21) 

t a =2 f - ?~r X dx 
Jo Sill <t> 

(22) 


For the dual-rotating propeller operating at the design con- 
ditions, the terms containing w are small, as is the case with 
the single-rotating propeller, and a close approximation to 
the drag loss is obtained if these terms are neglected. Fur- 
thermore, if it is assumed that the average of the resultant 
velocity W for the dual combination is equal to H 7 for the 
single propeller, equations (20) to (22) may be used for 
the dual-rotating propellers. Of course, for conditions other 
than the design condition, especially for very heavy loadings, 
exact drag-loss calculations require that the exact equations 
be used for either single-rotating or dual-rotating propellers. 

In summary, the equations for obtaining the propeller 
performance are given by the quantities c„ e, and P c and 
the drag-loss factors are given by t r and t a . 

The net thrust power is 


c« r — c« t a 

(23) 

The power input is 


Pc T = c»-f e "M r=Pc J T U 

(24) 

The efficiency is 


Cm la C ‘T 

^ P c+h P c T 

(25) 


where from equation (16) 

P c = 2 kw(1+w) 

The total power is also given by 

P=\ pFWP <r (26) 

It should be remembered that the . calculation is based on 
a given w. This procedure may seem unjustifiable since 
this parameter is not given by the specification but is the 
end result of a calculation based on the original data. The 
induced loss does not depend on the total-power coefficient 
P Ct but actually depends only on P c , and the quantity w 
cannot be obtained from the total-power coefficient. How- 
ever, the value of P Cr in most cases exceeds P c by not more 
than 2 percent or 

P e =0.98 P CT 


Since P c in equation (16) is based on w and the diameter 
of the final wake and since the value of P tF in equation (24) 
is based on the propeller diameter which is slightly larger 
than the diameter of the final wake, a veiy close approxima- 
tion to W is usually given by equation (16). Therefore, 

Pc T ~Pc= 2kw(1 4-w) ^1 

In some cases it may be necessary to calculate t, to obtain n 
more exact value of P C) especially if the blade profile drag is 
large. 

PROCEDURE FOR DESIGN OF PROPELLER 

FIGURES USED IN PROPELLER DESIGN 

The information necessary to design a propeller for any 
operating condition is given in the figures. Figure 3 gives 
the circulation function K{x) interpolated for even fractions 
for 2-, 3-, 4-, 6-, and 8-blade single-rotating propellers. The 
circulation function for the 2-blade propeller was taken 
directly from reference 5; for the 3-blade propeller, from 
reference 6; and for the propellers having a greater number of 
blades was recalculated from the Goldstein tip correction 
factors as given in reference 7. Figure 4 gives K(x) for dual- 
rotating propellers with 4, 8, and 12 blades. These values 
for the dual-rotating propellers were taken from data of ref- 
erence 1. Figure 5 gives the mass coefficient « for various 
numbers of blades for single-rotating propellers. Figure 0, 
which was taken from reference 1, gives k for dual-rotating 
, propellers. The ideal efficiency vt is plotted against ti; for n 
range of values e/ir iu figure 7, against c,/k in figure 8, and 
against PJk in figure 9. The data for figures 7 and 8 were 
taken directly from reference 4 and the data for figure 9 were 
recalculated by the use of equation (16) and figure 7. Fig- 
ures 7 to 9 apply to either single- or dual-rotating propellers. 
The propeller efficiency may be calculated from either of these 
figures ; however, in this report the efficiency has been deter- 
mined from P c /k as given in figure 9. 



V+Ifi 

Figure 5,— Mass coefficient * against tor various numbers of blades for single-rotating 

propellers. 
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Fiocse 6. — Measured values erf mass coefficient * for dual-rotating propellers with various 
numbers of blades. Two-blade single-rotating propeller Included for comparison (refer- 
ence 1; 


Figure 10 gives values of e, k, and e[* : for 2- and 4-blade 
single-rotating propellers and figure 11 gives values for 4-, 
8-, and 12-blade dual-rotating propellers. The values of « 
for a propeller vrith a finite number of blades have not pre- 
viously been published, but the values of e and e/g for an in- 
finite number of blades are given in figure 4 of reference 4. 
The method for calculating e/g and e is given in the following 
section. 

PROPELLER SELECTION 

In the selection of a propeller for a given airplane installa- 
tion, the engine power, the forward speed, and the design 
altitude are usually specified. The selection consists of the 
determination of the number of blades, the propeller solidity, 
the propeller diameter, and the rotative speed. The ideal 
propeller efficiency for any combinations of these variables 
can be readily obtained with the use of the charts. The 
procedure for a given blade number, propeller diameter, and 
rotative speed for either single or dual rotation is as follows: 
First, calculate the totaI-pow r er coefficient 

P = 

r C T 1 — 



Figcbk 7.— Propeller efficiency against w (reference 4). 


and then use this value for the ideal coefficient 

Pcj , « P c = 2 kw (1 + w) ^ 

to find w. 


It was shown in reference 4 that the dependence of the 
efficiency on e/s in the efficiency formulas is very small and 
that it is sufficient to know only the approximate value of 
e/s. An exa min ation of the formulas for e, and P c shows that 
their dependence on e/s is also small. It was further con- 
cluded in reference 4 that e/s is only slightly greater than s 
and that the practice of using e/s instead of s is considered 
satisfactory for design purposes. However, there appears 
to exist a simple relation between the axial-loss factor e and 
the total-loss factor s. This relation takes on the form of a 
differential equation 


e , 1_ X ds 

s“ 1+ 2sdX 

where 

% 1 V+w 

X 7lD 0 

This relation has been checked and found to be exact for 
an infinite number of blades, and numerical checks for a 
two-blade propeller were in very close agreement. It is 
considered accurate for an empirical relation for design 
purposes for propellers of other numbers of blades. 
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c,/k 

FioraE 8.— Propeller efficiency against c,[k (reference 4). 


First obtain 

x -IX. 

as a first approximation to X for use in the calculations. 
Then read off k and d<jd\ from the appropriate charts of k 

V _ 

against ^ (1-f-w) for several values of w (figs. 5 and 6). 


Curves of e, k, and e//c are plotted against -jy (l-J-tS) in fig- 
ures 10 and 11. Next plot a curve for the right side of the 
equation for P c against w. Where tliis curvo intersects the 
horizontal line, P c —Pc T is the desired point. This value may 
be checked from the chart by inserting the values obtained 
from the plot in the equation. Thus are obtained k, w, 
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^ (l+i/5), and. e/ic. From the chart of PJk (fig. 9), the 

optimum efficiency may be obtained. 

The following examples illustrate the method of deter- 
mining the optimum distribution of bc t along the radius for 
both single-rotating and dual-rotating propellers that give 
the ma ximum possible efficiency (neglecting blade profile 
drag) that can be obtained with either propeller for one 
specified design condition. 

ELLUSTEATIVE EXAMPLES 

Single rotation. — Let the following data specify the 
propeller design conditions: 


Power, horsepower 2,000 

Density, slugs per cubic foot 0.001065 

Velocity, miles per hour 425 


The propeller selection has been made to the extent that the 


following data specify the propeller: 

Rotational speed, n, revolutions per second . 23 

Propeller diameter, D, feet 12 

Number of blades, B 4 

VI nD 2.258 


The total P CT from the given conditions is 





(2000) (550) 

1(0.001065) (623)M6} ! 
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HH 

iim 'iiiiflfcn iniiiil 


~nD 

0 * ^PUT <b) rour-WadeslaRK-iotatiWProtvlVf. 

165al p,_<,.075. 

1 A plot of P« against w gives a value of ” 

The value of P, D and should be used I Then, y (2.258) ( 1 .155) =2.01 

instead of on^e propeher^a ^ may be ob- iZ5 (1+W) ^ 

to calculate w. approximations but the two e e y __ _2 gi and the 

»»*- 10 *»•■*"* ** 

® “ d Pt * gi ™ b7 ssvssft " *- 9 - 

equation (16) as *=0.201 


P e = 2kw( 1+^) (l+f W ) 


where 


P e =P e 


B ™lues of « are selected to «**£ "ft* 2TS 
Sot^"Sd e <or the four-blade emgle-rotattn* 
propeller : -r I 


-r,ss;s=s-- ~ 

k = 0.201 

*J= 0.373 

K 

and ,,,=0.929 

1 . , fbe single-rotating propeller 

Z£%3£'**** " om *“ a,iou co) 

, . sin*</> _ 

1+v — = 2u’K(x) 

■ ’ c, ”(Tl|®)(i+P ic0B ’ 4 ') 
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nt) tlD 

(a) Four-blade dual-rotating propeller. Cb) Eight-blade dual-rotating propeller. 

Floras II.— Values of > and c for dual-rotating propellers. 


Values of the circulation function K(x) at each station are 

y 

obtained from figure 3 (b) at (l+w)=2.61 and the angle 


of the relative velocity at the propeller is given for each 
station by 


tan 4>= 


1 V 
x tlD 


X 


Performing these calculations for w= 0.155 gives the values of 
<rC[ and bc t in the following table (the blade-width distribu- 
tion, in feet, for a constant c t of 0.5 is also given) : 


X 

tan 4 

RW 

•■Cr 

bet 

b 

m 

bet 

(bet) I.JM 

ai 

7.74 

0.033 

0.0S42 

0.079 

0.158 

a 167 

.2 

3.870 

. 078 

.0067 

.182 

.364 

.386 

.3 

a. 580 

.133 

.1054 

.298 

.596 

.631 

.4 

1.935 

.185 

.1044 

.393 

.786 

.833 

.fi 

1.543 

.225 

.0952 

.449 

-80S 

.0 2 

. S 

1.200 

.260 

.0855 

.483 

.966 

1.023 

.7 

Lioe 

.271 

.0716 

.472 

.944 

L00Q 

.8 

.968 

. 257 

.0654 

.417 

.834 

.880 

.9 

.860 

.204 

.0364 

.300 

.SIS 

.655 

.95 

.815 

.146 

.0241 

.216 

.432 

.453 


Dual rotation. — The procedure is repeated for a 12-foot- 
diameter four-blade dual-rotating propeller for the same 
design conditions as used for the single-rotating propeller. 
The following table is obtained for the four-blade dual- 
rotating propeller ^values of k and e//c were found from 
figure 11 (a)): 


w 

K 

«/* 

p. 

0 

a 472 

a 580 

0 

.1 

.432 

.547 

.1002 

.2 

1 .398 

.619 

.211 


In this case a plot of P e against w gives a value of u?=0.075 
at P c = 0.075. Therefore, 


and 


(l+w) =2.426 
*=0.442 


^=0.170 

K 


ij<=0.964 


024778- — 51 8 
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(e) Twelvo-blade dual-rotating propeller. 


Fiourb 11.— Concluded. 


It is seen that the important parameter, the mass-flow 
coefficient, is 0.442 for the dual-rotating propeller and is only 
0.201 for the single-rotating propeller. The efficiency 
(without drag) is 96.4 percent for the dual-rotating propeller 
but is only 92.9 percent for the single-rotating propeller. 

For the dual-rotating propeller the values of <rc t may be 
found for the front component from equation (12); thus, 


(<tCi)f — 


V _1_ {l -\-w)w sin (ftp 
nD 

1 +£ kid sin a 0a 


K(x) 


and for the rear propeller from equation (13) 

, .. n V 1 (l+«j)«i sin (ft, TJ r,^ 

((rc ‘ )s ~w « -.3 “ K(x) 

1 kiv sin^o 

Equation (9) gives <f> F by 

tan ^ [l +1 v (l +1 k tanv)] 


and <t>B is given in equatiou (11) by 

tan ^ [l ■ (l ■ k tanv)] 

Values of the circulation function K(x) are obtained from 

V-\-w V 

figure 4 (a) at the appropriate value of — p (l+t5). 

Performing these calculations for V=623, n=23, V— 12, 
w= 0.075, and k= 0.442 gives the values of tan <j> and <rc { in 
the following table: 


X 

ZW) 

tan 4>r 

tan** 

(<rci)r 

frci)* 

(ba)r 

<6f0* 

0.1 

0.620 

10.708 

4.148 

0.3510 

0.3460 

0.064 

0.084 

.3 

.807 

2.608 

2303 


.1058 

.524 

.618 

.4 

. 696 

L 916 

1.812 

.0743 

.0738 

.864 

.550 

.6 

.672 

L 618 

1.465 

.0532 

.0628 


.190 

.8 

.838 

1.258 

1.227 


.0369 

.422 

,418 

.7 

.488 

1.075 

L 056 


.0285 

.378 

. 370 

.8 

.417 

.639 

.920 


.0199 

.302 

.301 

.9 

.317 

.833 

.824 

.0126 

.0126 

.214 

.214 

.96 

.241 

.789 

.781 

.0038 

.0088 

.167 

.157 


A comparison of the optimum distribution of bc t along the. 
blade for the dual-rotating propeller from this table with, 
the optimum distribution for the single-rotating propeller 
as given in the preceding section shows that, if approxi- 
mately constant c, is absorbed along the blade, wide dif-. 
ferences in blade plan form will result for the two propellers- 
designed for the Bame operating condition. For the operat- 
ing conditions selected, the maximum bc t for the single-: 
rotating propeller occurs near the 0.6 radius and tapers rapidly 
toward the tip and the hub, being only slightly over 10. 
percent of its maximum value at the 0.1 radius. On the. 
other hand, the minimum value of bc t for optimum distribu- 
tion for the dual-rotating propeller occurs at the propeller 
tip and progressively increases toward the inner radii. The 
value of bci at the 0.1 radius is four times its value at the 
0.96 radius. 

Since the design of the dual-rotating propeller calls for 
high loading over the inner sections, the efficiency of the dual- 
rotating propeller is less susceptible to compressibility losses 
which normally occur near the propeller tip for operation at 
high tip Mach numbers. The compressibility losses may 
be reduced by reducing the width of thoso sections or by 
reducing the operating lift coefficient. 

Effect of blade drag on efficiency. — Tho loss in efficiency 
due to the profile drag of the blades can bo calculated from 
equations (20) to (22) if the blade-width distribution and 
profile-drag coefficients at the operating c, arc known. Inas- 
much as structural requirements may determine the shape 
of the blade, especially over the inner radii, only ono example 
is given. The equations, however, may be applied to any 
plan form. The example selected is for the four-blade 
single-rotating propeller on which the induced efficiency has 
beeDL previously calculated. The shank sections of tho 
propeller blade were assumed to bo round, similar to thq 
Hamilton Standard Propeller No. 3165-6 and the blado 
plan form from r=0.3 to *=1.0 was mado optimum for a c\ 
of 0.5. The profile-drag coefficients for the several radii are 
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the same as given in reference 7 for the Hamilton Standard 
Propeller No. 3155-6 which has Clark Y sections and are 
given in the following table. It is assumed that a spinner 
covers the inner 0.2 of the radius. The distribution of 
ac t with x and of sin <j> with x have been included in the table: 


X 

<rci 

r 

Cd 

sintf 

—X 

sin* 

sin 

0.2 

0.0967 

a 1934 

0.400 

0.96S 

a 01600 

0.00064 

.3 

. 1054 

.2108 

.100 

.032 

.00697 

.00061 

.4 

.1044 

.20SS 

.020 

.889 

.00188 

.00030 

.5 

.0952 

.1904 

.010 

-W0 

.00113 

.00028 

.6 

.0855 

.1710 

.008 

.790 

.00104 

.00037 

.7 

.one 

.1432 

.007 

.742 

.00093 

.00046 

.8 

. 0554 

.1106 

.006 

.696 

.00077 

.00049 

.9 

.0364 

.0728 

.006 

.652 

.00060 

.00049 


Performing the integrations and substituting in the for- 
mulas gives for rotational-drag-loss coefficient 


o ris> 

<r-r=i ; 

A* Jo . 2 ! 


UCj 

0.2 sin <ji 


i 8 dx 


= oil 6 (°- 000348 )= 0 - 0014 

and for the axial-drag-loss coefficient 

U=2 r^xdx 

J o.2 sm 4> 

=2(0.00213) =0.0043 

The induced thrust coefficient has been given by equa- 
tion (14) as 

c,=2tiw |^1 +w Q+f)J 

=2(0.201) (0.155) ^1+0.155 (J + 0.29^=0.0700 
and the induced power coefficient by equation (16) as 
P c =2kw(1+w)(\ +~ 

=2(0.201) (0.155) (1.155) (1.045) =0.0754 
The induced efficiency is 

e. 




0.0700 „ 
- 03754- 0 - 929 

With drag included, the total thrust is given by 
= 0.0700 -0.0043 = 0.0657 


and 

Pc T =P'+tr 

=0.0754+0.0014=0.0768 


The efficiency is 



0.0657 

0.0768 


=0.855 


Thus it is seen that the blade drag of the magnitude given 
in the preceding table reduces the propeller efficiency from 
92.9 percent to 85.5 percent for the propeller operating 
conditions given. 


CONCLUDING REMARKS 

A comparison of Theodorsen's propeller theory with the 
conventional vortex theory shows that the optimum load 
distribution along the blade for single-rotating propellers 
obtained by the two theories is essentially identical and as a 
result the op timum efficiencies are the same for a given 
operating condition. Theodorsen's theory has the advantage, 
however, that the optimum efficiency for any design con- 
dition can be obtained quickly and accurately by the use of 
the mass coefficient k without any laborious calculations and 
before the final design is made. 

The distribution of the circulation function K(x) for the 
idealized dual-rotating propeller is radically different from 
the existing values for the single-rotating propeller that 
have been previously used for the dual-rotating propeller. 
Also, the mass coefficient ic for the dual-rotating propeller 
is larger than the sum of the values for two single-rotating 
propellers. These quantities, which are not available from 
mathematical computations but are obtained from the 
electrical-analogy method of Theodorsen, are used herein 
for obtaining the opt imum load distribution along the blade 
for the dual-rotating propeller. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., March 15, 1948. 
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